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ABSTRACT 

Massive early-type galaxies conimonly have gas discs which are kinematically mis¬ 
aligned with the stellar component. These discs feel a torque from the stars and the 
angular momentum vectors are expected to align quickly. We present results on the 
evolution of a misaligned gas disc in a cosmological simulation of a massive early-type 
galaxy from the Feedback In Realistic Environments project. This galaxy experiences 
a merger which, together with a strong galactic wind, removes most of the original 
gas disc. The galaxy subsequently reforms a gas disc through accretion of cold gas, 
but it is initially 120° misaligned with the stellar rotation axis. This misalignment 
persists for about 2 Gyr before the gas-star misalignment angle drops below 20°. The 
time it takes for the gaseous and stellar components to align is much longer than 
previously thought, because the gas disc is accreting a significant amount of mass for 
about 1.5 Gyr after the merger, during which the angular momentum change induced 
by accreted gas dominates over that induced by stellar torques. Once the gas accretion 
rate has decreased sufficiently, the gas disc decouples from the surrounding halo gas 
and realigns with the stellar component in about 6 dynamical times. During the late 
evolution of the misaligned gas disc, the centre aligns faster than the outskirts, result¬ 
ing in a warped disc. We discuss the observational consequences of the long survival 
of our misaligned gas disc and how our results can be used to calibrate merger rate 
estimates from observed gas misalignments. 

Key words: galaxies: evolution - galaxies: formation - galaxies: elliptical and lentic¬ 
ular, cD - galaxies: kinematics and dynamics - methods: numerical 


1 INTRODUCTION 

Elliptical and lenticular galaxies are thought to form at high 
redshift throu gh repeated mergers of smaller galaxies (e.g. 
lToomrelll977l l. Such ‘early-type galaxies’ (ETGs) are gener¬ 
ally among the most massive galaxies in our universe. Lo¬ 
cal ETGs are mostly gas free and dominated by old stars, 
while lower mass spiral galaxies have a more massive inter¬ 
stellar medium (ISM) and are still forming stars. This anti- 
hierarchical growth of luminous structure (termed ‘down- 
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sizing’; iGowie et al.l Il996ll is one of the key puzzles in an 
otherwise hierarchical cold dark matter (CDM) universe. 


The gas-poor nature of ETGs makes them excellent lab¬ 
oratories in which to study the effects that stellar mass- 
loss, cooling from the hot halo, gal axy mergers, a nd accre¬ 
tion have on galaxy evolution (e.g. lKaviraill2014h . Studies 
of statistically complete samples of ETGs have shown that 
about a quarter of these objects have more t han 10^ M(r, 
of co l d molecular gas dCombes et Tl l2007l : IWelch et al'l 
I 2 OI 0 I : lYoung et akl 1201 il l and appro ximately 40 per cent 
have sizeable atom i c gas reservoirs llMorganti et al.l l2006l : 
I Sage fc Welchllioo^ : IsAa et'aDl2012l~). Dust is also p resent 


in about 50 per cent of ETGs (e.g. Smith et al.ll2012 1. This 
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ISM mat erial appears to be regeneratiiiK previously gas-poor 
galaxies J Young et alJ 12014 1. but its origin is debated. 

Observationally, in field environments the gas and dust 
in ETGs seems to come mainly from ex ternal sources, such 
as mergers and c o ld gas accretion (e.g. Kna£ 2 ._et_^ 198^ 
Sarziet_aJj_ 200d McDermid_et_^ 200d; Oosterloo^£_^ 


2010 l : lDavis et al1l201ll : ISinith et aLll2012l : IDuc et aLll2015ll . 

Theoretically, however, it has been argued th at the main 
mode of gas supply is cooling from hot halo gas (jLagos et al.l 
[ 20 T 3 , [mi), with minor mergers playing a secondary role. 
These semi-analytic models had to make assumptions about 
some of the relevant physical processes, which need to be 
verified by hydrodynamic simulations. 

Here we concentrate on understanding the kinematic 
decoupling which has been observed bet ween the gaseous 
and stellar components of ETGs (e .g. ISarzi et al.l l 200 fil : 
IPavis et al.ll 201 ll : iKatkov et al.ll20l4 l. About 50 per cent 
of local field ETGs have gas kinematically mis a ligned from 
their stellar kinematic axis ([Pavis et al.l I 2 OI 1 I : ISerra et al.l 
120141 . There are indicat ions of the same behav iour at high 
redshift. At 2 « 1 — 2.5, IWisnioski et al.l ll2015l l found large 
misalignments between the gas kinematic and stellar pho¬ 
tometric axis in several galaxies. These galaxies tended to 
have redder colours and lower specific star formation rates 
than the majority of galaxies in their sample, consistent with 
these objects being ETGfEI 

A priori, it seems likely that that the origin of mis¬ 
aligned material has to be external, as gas produced by 
internal processes (e.g. stellar m ass-loss) should sh are the 
angular momentum of the stars llPavis et al.l 1201 il l . Thus 
the observed incidence of misalignments gives a lower limit 
for the fraction of galaxies where the gas must have been 
supplied by external processes. However, in an axisymmetric 
potential the gas will over time relax back into one of the pre¬ 
ferred axes, becoming exactly co- or counter-rotating. This 
is because the gas feels a torque in the (quadrupolar) grav¬ 
itational potential of the stars. This torque will make the 
gas disc precess around the angular momentum direction of 
the stars. Because the torque is larger at smaller radii, gas 
will precess at a faster rate. Pue to the fact that the gas is 
not on exactly circular orbits, there will be collisions of gas 
clouds precessing at different rates. These collisions work to 
realign the differentially precessing gas ‘rings’. Pissipative 
forces destroy the component of the angular momentum of 
the gas perpendicular to the stellar angular momentum di¬ 
rection, until they are aligned or 180° misaligned, the most 
stable configurations. 

Theoretical studies of the relaxation of misaligned gas 
discs in the potential of elliptical galaxies have considered 
this effect in both axisymmetric and triaxial halos and found 
that in both cases the a lignment process typically takes a few 
dynamical times Ifriy,,: iTohline et al.lll982l : iLake fc NormarJ 
ll984i . lLake fc NormanI 1 1983ll found the relaxation time, i.e. 
the time it takes a misaligned disc to settle into the plane, to 
be similar to the differential precession time, which is equal 
to treiax ~ tdyn/e, where £ is the eccentricity of the potential. 


For a typical lenticular galaxy e « 0.2 (iMendez-Abreu et al.l 
l2008ll . and Geiax ~ 5tdyn. This time is very short in the 
centre of local ETGs (about 100 Myr), and thus very few 
misaligned objects should be observed. It is this puzzling 
mismatch between the theoretical predictions of fast relax¬ 
ation and observations of many misaligned discs that we aim 
to investigate here. Pavis & Bureau (in preparation) point 
out that if one can understand the time this relaxation pro¬ 
cess takes, one can obtain (at low redshift, where mergers 
are expected to dominate the supply of cold gas) an estimate 
of the rate of gas-rich mergers in massive galaxies, a quan¬ 
tity which is important to galaxy evolution in a standard 
AGPM framework. 

In this paper we conduct a case study, investigating the 
evolution of a gas disc embedded in the centre of a large 
ETG within a cosmological, hydrodynamic simulation. This 
gas disc is accreted in a misaligned configuration and relaxes 
slowly back towards the stellar kinematic axis as the galaxy 
evolves. We find that the gas-star misalignment persists for 
more than 2 Gyr, much longer than previously expected, 
due to continuous gas accretion. In Section [5] we describe 
the simulation we used. In Section |3] we present our results, 
with Section 13. 1 1 focussing on the galaxy’s accretion history 
and metallicity and Section 13.21 focussing on the gas disc’s 
misalignment and relaxation. The main results are shown 
in Figures [S] and [7] We discuss our results and conclude in 
Section U) 


2 METHOD 

This work is par t of the ‘feedback in realistic environments’ 
(FIRE) projecl|3 llHopkins et al.l[20l4 l. which consists of sev¬ 
eral cosmological ‘zoom-in’ simulations. Here we make use of 
one galaxy, ml3, with stellar mass Mstar = 10^^ Mq, simi¬ 
lar to that of giant ETGs at 2 ; = 0, and with a halo mas^ 
A/haio = 10*^® Mp) . Its initial cond itions are taken from the 


AGORA project (IKim et alJ 

2014 

). The simulation is fully 

described in iHonkins et al. 1 

2014 

) and references therein. 


but we will summarize its main properties here. 

We use a version of TreeSPH (P-SPH) which adopts the 
Lagrangian ‘pressure-entropy’ formulation of the smoothe d 
particle hydrodynamics (SPH) equations (iHopkinsI l2013l l. 
The gravity solver is a heavily m odified version of GADGET-2 
(last described in ISpring^l2005I) . P-SPH also includes sub¬ 
stantial improvements in the artificial viscosity, entropy dif¬ 
fusion, adaptive timesteppin g, smoothing kernel, and gravi¬ 
tational softening algorithm llHopkins et al.ll20l4l . 

A ACDM cosmology is assumed with parameters consis¬ 
tent with the 9-yr Wilkinson Microwav e Anisotropy Probe 
(WMAP) results llHinshaw et al.|[2013h . Hm = 1 — Ha = 
0.272, Hb = 0.0455, h = 0.702, as = 0.807 and n = 0.961. 
The (initial) particle masses for dark matter and baryons are 
2.3 X 10® M© and 4.5 x 10® M©, respectively, for our fiducial 
resolution. The minimum physical baryonic force softening 
length for our fiducial resolution is 28 h~^pc. We adopt a 
quintic spline kernel with an adaptive size, keeping the mass 


^ Star-forming galaxies have much higher gas fractions than 
ETGs and are not likely to exhibit significantly misaligned gas 
discs, because they are unlikely to accrete enough material (and 
angular momentum) to tilt the entire gas disc. 


^ http://fire.northwestern.edu/ 

® We define halo mass as the total mass inside a spherical region 
within which the mean density is 200 times the mean density of 
the Universe 
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within the kernel approximately equal. The average number 
of neighbours in the smoothing kernel is 62. The contamina¬ 
tion fraction, i.e. the fraction of the total dark matter mass 
that is in lower resolution (higher mass) particles, at 2 = 0 is 
0.07 per cent within 20 kpc, the region we study. Therefore, 
we do not expect contamination to affect our results. 

Star formation takes place in molecular, self-gravitating 
gas above a hydrogen number density of rig > 10 cm“®, 
where the molecul a r fra ction is calculated following 
iKrumholz fc GnedinI l^2011^ and t he self-gravitating crite¬ 
rion following Hopkins et al.l ll2013l l. Stars are formed at the 
rate p* = Pmoiecuiar/tff, where tg is the free-fall time. Star 
particles inherit their mass and metal abundances from their 
progenitor ga s particle. We assume an initial stellar mass 
function from iKroupal (l2002l 'l. Radiative cooling and heating 
are computed in the presence of the cosmic microwave back¬ 
ground (CM B) radiation and the ultrav iolet/X-ray back¬ 
ground from iFaucher-Giguere et aP (l2009l 'l. Self-shielding is 
accounted for with a local Sobolev/Jeans length approxima¬ 
tion. We impose a temperature floor of 10 K or the CMB 
temperature. 

The primordial abundances are X = 0.76 and Y = 
0.24, where X and Y are the mass fractions of hydrogen 
and helium, respectively. The simulation has a metallic- 
ity floor at Z = 10“"^ Zq, because yields are very uncer¬ 
tain at lower metallicities. The abundances of 11 elements 
(H, He, C, N, O, Ne, Mg, Si, S, Ca and Fe) produced by 
ma ssive and intermediate - mass stars are computed f ollow¬ 
ing Iwaniptoet_alJ lll999l l , IWooslev fc WeavCT lll995l ~l , and 
lizzard et al.l ( 20ol i. Mass ejected by a star particle through 
stellar winds and Type la and Type H supernovae is trans¬ 
ferred to the gas particles in its smoothing kernel. 

The FIRE simulations include an implementation of 
stellar feedback by supernovae, radiation pressure, stellar 
winds, and photo-io nization and photo-electric heating (see 
iHopkins et al.ll2014l and references therein for details). For 
the purposes of the present paper, we emphasize that these 
simulations produce galaxies with stellar masses and metal¬ 
licities reasonably consistent with obs ervations over a wide 
range of dark matter halo masses llHopkins et al.l I20l4 
iMa et al] l2015h . This is a conseque nce of galactic wind s 
driven by stellar feedback (see also iMuratov et al.l 1201511 . 
Feedback from active galactic nuclei is not included. 

All distances in this work are given in proper kpc. Many 
results are shown as a function of lookback time, tiookback. In 
our cosmology, tiookback = 0, 1, 2, 3, 4, and 5 Gyr correspond 
to 2 = 0, 0.075, 0.16, 0.26, 0.37, and 0.49, respectively. A 
major galaxy merger occurs at tiookback = 3.4 Gyr (2 = 0.3), 
which is indicated in the figures by the grey shaded region. 


3 RESULTS 

In the present work, we study the evolution of the gas 
disc embedded in the centre of a large early-type galaxy. 
Specifically, the galaxy experiences a gas-rich merger at 
fiookback = 3.4 Gyr (2 = 0.3), which pushes out most of its 
gas disc. The original disc material reaccretes quickly, but 
is soon thereafter permanently removed from the system in 
a violent outflow. When the galaxy subsequently regrows 
a gas disc, this disc is initially misaligned with the stellar 
rotation, but eventually relaxes into the plane of the stars. 


To illustrate this evolution, we show the average line- 
of-sight velocity of all the gas in a (20 kpc)® box around the 
galaxy centre in Figure [T] at several lookback times as indi¬ 
cated in each panel. The colour scale runs from —300 km s“^ 
(blue) to 300 km s“^ (red) and is saturated at the extremes. 
The image is rotated such that the stellar ‘disc’ is edge on 
and horizontal, i.e. its angular momentum vector points in 
the j/-direction. The shape of the stellar component is shown 
as white contours which connect areas of equal mass surface 
density, with logarithmically spaced contour levels. This also 
shows that the galaxy is an ETG during the entire evolu¬ 
tion we consider here. At tiookback ^ 4 Gyr, the stellar axis 
ratio is reasonably small, resembling an SO. At later times, 
after the merger, the axis ratio is larger and the galaxy more 
resembles an elliptical. 

The Erst panel shows the gas disc 0.7 Gyr before the 
merger, the second panel right after the merger when most 
of the gas disc has been removed. The galaxy reforms a gas 
disc through accretion, but it is misaligned with the stars, 
as can be seen in the third panel. The fourth panel shows 
that the disc tilts and precesses with respect to the stars, 
which is due to continued gas accretion as explained below. 
The central kpc precesses faster than the outskirts, starting 
a warp. At late times, it develops a warp in the outer parts 
as well, as shown in the fifth panel. Here, the evolution is 
driven by stellar torques, which are stronger in the central 
part than in the outer parts, resulting in the centre aligning 
first. By the end of the simulation, shown in the last panel, 
the gas disc is almost completely aligned with the stars. 


3.1 Evolution of galaxy properties 


At tiookback = 5 Gyr, the galaxy’s stellar component shows a 
small disc structure in its kinematics, although it is clearly 
bulge-dominated. The stellar disc component disappears af¬ 
ter the last major merger, but there is still a preferred ro¬ 
tation direction. The specific angular momentum decreases 
after the merger. We measure the mass-weighted 2D edge-on 
axis ratio, b/a, and Xr, an observable proxy of the specific 
angular momentum dEmsellem et al.ll2007l . l201l| j. 


Xr 


R2d|ulos| 


( 1 ) 


where R 2 D is the half-mass radius measured in projection, 
ulos is the line-of-sight velocity, and a is the velocity disper¬ 
sion in each (0.1 kpc)^ pixel. We find that this galaxy is a fast 
rotator both b efore and after the me rger, according to the 
classification of lEmsellem et akl ll201lll . Xr > 0.31^1 — b/a. 
Before the merger, Xr = 0.4—0.5. Afterwards, Xr decreases, 
mostly due to the reaccretion of the stars removed in the 
merger, until it levels off at about 0.25 at fiookback < 2.6 Gyr. 
It is possible that the luminosity-weighted Xr would be even 
higher, since young stars have higher velocities, so the clas¬ 
sification of this galaxy as a fast rotator is robust. The prop¬ 
erties of the gas disc match observations of local ETGs well 
(see the discussion in Section |4]). 

Figure [ 2 ] shows the total gas mass (solid, black curve) 
and the star-forming gas mass (dashed, blue curve), i.e. the 
gas at riH > 10 cm“®, in the central 5 kpc and the stellar 
mass (dotted, red curve) in the central 20 kpc. The stellar 
mass inside 5 kpc is only 0.2 dex lower than the total stellar 
mass shown. The mass of this galaxy is dominated by its 
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Figure 1. The average line-of-sight velocity of all gas particles in a (20 kpc)® box centred on the galaxy at several lookback times, i, 
as indicated in each panel. The colour scale runs from —300 (blue) to 300 km s“^ (red) and is saturated at either end. The image is 
rotated such that the stellar angular momentum vector points in the j/-direction (i.e. upwards). The shape of the stellar component is 
shown as white contours which connect areas of equal mass surface density, with logarithmically spaced contour levels. The first panel 
shows the gas disc before the merger, the second panel right after the merger when most of the gas disc has been destroyed. The galaxy 
reforms a misaligned gas disc through accretion, as can be seen in the third panel. This gas disc is seen to rotate in the fourth panel, 
due to continued misaligned gas accretion. At late times it exhibits a warp, clearly visible in the fifth panel, because the gas disc torques 
with the stars and the stellar torques are stronger in the centre, which therefore realigns first. The last panel shows that the gas disc is 
almost completely aligned with the stars. The gas disc remains misaligned for much longer than a few times tdyn- 


stars within 25 kpc, outside of which dark matter dominates. 
At a = 0, the stellar mass inside 20 and 5 kpc is 10^^ and 
6 X 10^° Mq, respectively, and has only grown by a factor 
of 1.7 and 1.3 in the past 5 Gyr. In the central 5 kpc, about 
half the stellar mass increase is due to the accretion of stars 
originally belonging to the merging galaxy and the other 
half is due to star formation. The stellar half-mass radius at 
late time is 3.4 kpc for stars within 20 kpc. 

The galaxy that merges and destroys most of the gas 
disc does so on its second passage. Just before the merger, 
the satellite’s gas mass inside 5 kpc and stellar mass inside 
20 kpc are 2.2 x 10® and 2.2 x 10^® Mq, respectively. This 
is a factor of 1.4 larger than the central galaxy’s gas mass 
and a factor of 3.8 smaller than the central galaxy’s stellar 
mass at the same time. With a stellar mass ratio of 1:4 this 
merger is on the border of being classified as a major merger, 
whereas in gas mass the merger ratio is close to 1:1. 

The gas mass decreases by a factor of a few after the 
merger to 5 x 10® Mq. Most of the gas in the outer parts 
of the disc is removed, but the central kpc remains rela¬ 
tively intact. The half-mass radius of gas within 5 kpc is 
only 0.4 kpc right after the merger, but grows steeply af¬ 
ter fiookback = 3 Gyr, staying between 2.5 and 3.2 kpc after 
tiookback = 2.1 Gyr. The total gas mass continues to grow 


down to fiookback = 1-7 Gyr to a total of 3 x 10® Mq and 
stays stable after that. The star-forming gas mass also de¬ 
creases after the merger, but does not build back up as much 
as the total gas mass. Additionally, it decreases steadily af¬ 
ter fiookback = 1-7 Gyr, due to star formation and relatively 
low accretion rates. The star-forming gas mass is an order 
of magnitude lower than the total gas mass at the end of 
the simulation, whereas at e.g. fiookback = 4 Gyr the differ¬ 
ence was less than a factor of three. The half-mass radius 
of the star-forming gas varies between 0.2 and 0.9 kpc over 
the past 5 Gyr. These gas masses and radii are consistent 
with the range observe d in ETGs (e.g. lYoung et al.l I2OI1I : 
iDavis et al.ll2C)lll. l2013l) . 

Inflow and outflow rates in the past 5 Gyr are shown 
in Figure [3] as the solid, black curve and dashed, light-blue 
curve, respectively. Additionally, the stellar mass-loss rate 
and star formation rate (SFR) are shown by the dot-dashed, 
purple curve and the dotted, red curve, respectively. All 
rates are measured within 5 kpc of the centre and averaged 
over the time between snapshots (which varies from 36 to 
69 Myr). We checked that there is no significant difference 
when using a SFR averaged over 100 Myr instead. The stel¬ 
lar mass-loss rate shown here is defined as the increase in gas 
mass in the centre not due to particles leaving or entering 
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Lookback time [Gyr] 


Figure 2. The total (solid, black curve) and star-forming 
(dashed, blue curve) gas mass within 5 kpc and stellar mass (dot¬ 
ted, red curve) within 20 kpc as a function of lookback time. The 
stellar mass is close to 10^^ Mq and only grows by a factor of 1.4 
over the past 5 Gyr. The total gas mass is a few times 10^ Mq for 
most of its recent history, but drops to 5 x 10^ Mq shortly after 
the merger (indicated by the grey band). The star-forming (i.e. 
dense) gas mass is highest before the merger (above 10® Mq) and 
decreases to 2 x 10® Mq after the merger. At late times, the total 
gas reservoir remains constant, but the star-forming gas reservoir 
decreases, due to star formation. 



Figure 3. The gas inflow rate (solid, black), gas outflow rate 
(cyan, dashed), star formation rate (dotted, red), and stellar 
mass-loss rate (purple, dot-dashed) within 5 kpc from the cen¬ 
tre as a function of lookback time. On average, the inflow rate 
dominates over the outflow rate and stellar mass-loss rate. After 
^lookback = 3.1 Gyr, there are two significant periods of gas ac¬ 
cretion, but there are no strong outflows. The gas accretion rate 
remains constant (and relatively low) at around 3 Mq yr~^ at 
^lookback ^ 1-3 Gyr. 

the inner 5 kpc. It is much lower than the gas accretion rate 
and thus subdominant at reforming the gas disc and fueling 
star formation after the merger. 

The galaxy experiences strong inflow and outflow be¬ 
tween tiookback = 4 and 5 Gyr, which is due to close flybys 
by satellites and to a minor merger. The rates are much 
lower afterwards, but peak at the time of the main merger 


at tiookback = 3.4 Gyr (at 48 and 59 Mq yr“^ for inflow 
and outflow, respectively) and shortly thereafter. After this 
time there are no major outflows. The galaxy experiences 
two episodes of relatively strong gas inflow (without associ¬ 
ated strong outflow), peaking at tiookback = 2.9 and 1.8 Gyr. 
The first peak is due to the accretion of the majority of the 
gas content of the galaxy that merges at tiookback = 3.4 Gyr. 
The gas accretion rate remains low in the last 1.5 Gyr, con¬ 
sistent with the lack of late-time mass growth as seen in 
Figure [2] 

The SFR of this galaxy at 2 = 0 is 1.4 Mq yr“^, which 
means its specific SFR is about 10~^^ y r~^, consistent wit h 
observed SFRs in gas-rich ETGs (e.g. IPavis et al JHoIi). 
This gal axy lies below the blue cl oud of star-forming galax¬ 
ies (e.g. iBrinchmann et ^ l2004l . but always has a non- 
negligible gas reservoir and residual star formation (see Fig¬ 
ures [2] and o. We also find a kinematically decoupled core 
in young sta rs, which are formed in the misaligned gas disc, 
as observed l|McDermid et af]l200fil l. The SFR fluctuates be¬ 
tween 0.1 and 7.0 Mq yr“^ over the last 5 Gyr. The average 
inflow and outflow rates over the past 3 Gyr, i.e. after the 
merger, are 3.8 and 1.7 Mq yr“^, respectively, and the av¬ 
erage SFR is 2.4 Mq yr“^. The average stellar mass-loss is 
only 0.6 Mq yr”*^ over the past 3 Gyr. 

Most of the gas that accretes onto the central 5 kpc in 
the past 5 Gyr does so in the ‘cold mode’, I.e. it never reached 
temperatures close to the vlrial temperature of the halo 
(Tvir ~ 1.4x 10® K) in its history. In fact, only approximately 
30 per cent of the gas that accretes after the merger ever 
had a temperature above 10®'® K, a temperat ure at which 


the coolin 

p time is relativelv short (see e.s. iKeres et al. 

2005 

2009 

: van de Voort et al. 201ll:lFaucher-Gieuere et al. 

2011 

van de Voort & Schavell2012^. The maioritv of the eas 


that reforms the disc therefore did not cool from a diffuse, 
hot, hyd rostatic halo, as wa s assumed in the semi-analytic 
model of iLagos et al.l (l2015ll . Instead, this gas was brought 
in with high density by the merging galaxy and via satellites 
that were stripped in the halo, i.e. the gas accretion is not 
smooth, but clumpy. 

A substantial fraction of the gas originally in the disc is 
removed during the merger, but reaccretes after only about 
0.1 Gyr. After another 0.1 Gyr a strong supernova-driven 
outflow occurs, reducing the gas mass once again, but the 
gas ejected in this event has not reaccreted by 2 = 0. Only 40 
per cent of the gas in the disc prior to the merger remains in 
the centre, forming stars. On the other hand, 2/3 of the gas 
brought in by the merging galaxy is accreted onto the central 
disc. Most of this accretion happens at tiookback = 2.9 Gyr, 
substantially increasing the total gas mass in the disc. 

The importance of gas accretion is also obvious when 
we look at the evolution of the average metallicity, [Fe/H], 
shown in Figure [T] for the gas within 5 kpc as a solid, 
black curve. [Fe/H] is the abundance ratio of iron to hy¬ 
drogen compared to that of the Sun and is defined as 
log^Q(A‘pe/VH)-logiQ(VFe/iVH)Q where Vpe and Ah are the 
number densities of iron and hydrogen. The stellar metallic¬ 
ity within 20 kpc and 3 kpc of the galaxy centre is shown 
by the dotted, red and dashed, purple curves, respectively. 
3 kpc is close to the stellar half-mass radius. 

The stellar metallicity within 20 kpc (3 kpc) rises 
slowly from [Fe/H]= —0.5 (—0.4) to —0.4 (—0.2) by 2 = 
0. The slightly subsolar stellar metallicity and the fact 
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Lookback time [Gyr] 


Figure 4. Mean gas metallicity within 5 kpc (solid, black) and 
stellar metallicity within 20 kpc (dotted, red) and within 3 kpc 
(dashed, purple) of the galaxy centre. The decrease in gas metal¬ 
licity after tiookback = 3 and 2 Gyr corresponds to periods 
of high gas inflow rates, as seen in Figure [3] This also shows 
that the galaxy is not reaccreting its pre-merger gas between 
^lookback = 3.2 and 1.7 Gyr, but instead is accreting new ma¬ 
terial with lower metallicity. At later times, the gas metallicity 
increases again, because the accretion rate is low enough that the 
enrichment by supernovae dominates. 



Lookbock time [Gyr] 


Figure 5. The misalignment angle between the angular momen¬ 
tum vector of the gas and stars as a function of lookback time 
is shown by the solid, black curve. A substantially misaligned 
disc is present for about 2 Gyr, much longer than a dynamical 
time. The misalignment angle for gas with radii between 5 and 
20 kpc, i.e. gas which will likely accrete onto the gas disc, is 
shown by the dotted, orange curve. The inner and outer misalign¬ 
ment angles evolve together after the merger, until they diverge 
3-t ilookback = 1.3 Gyr, after gas accretion slows down, enabling 
the stellar torques to realign the inner gas disc. 


that it decreases with increasing radius a re both consistent 
with the majority of observed ETGs (e.g. iKuntschner et al.l 
l2ni(t >. The gas metallicity decreases after the merger from 
[Fe/H]=0.6 to 0 by tiookback = 1-7 Gyr. This indicates that 
the accretion of fresh, metal-poor material dominates over 
the metals created by massive stars and expelled by stellar 
winds and is consistent with our finding that the gas mass 
increases substantially and gas accretion rates are high down 
to tiookback « 1.7 Gyr (see Figures [2] and [3|. Afterwards, the 
metallicity increases again to [Fe/H]=0.3 at « = 0, show¬ 
ing that at late times the metal dilution by accretion is less 
important than the metal enrichment by star formation. 

3.2 Misalignment and relaxation 

Figure [S] shows the misalignment angle between the angular 
momentum vector of the gas within 5 kpc and stars within 
10 kpc as a function of lookback time as a solid, black curv^f]. 
The misalignment angle for gas with i? = 5 — 20 kpc, i.e. 
gas which is likely to accrete onto the gas disc in the near 
future, is shown by the dotted, orange curve. The central 
disc shows no misalignment before the merger. The gas at 
7? = 5 — 20 kpc starts to misalign at tiookback = 4 Gyr, 
likely a result from gas stripping from close satellite flybys, 
but it does not affect the fairly massive central disc signifi¬ 
cantly. At the time of the merger, the disc loses most of its 
mass and subsequently accretes gas with misaligned angular 
momentum from the region around it. 

The gas at 77 < 5 kpc and 77 = 5 — 20 kpc show sim¬ 
ilar gas-star misalignment angles with the stars down to 

We use 10 kpc for stars instead of 20 kpc in order to be less 
affected by close satellite flybys, but the results are insensitive to 
this choice. 


tiookback = 1-3 Gyr. These angles decrease at the same rate 
and since the misalignment angle of the gas outside the disc 
is slightly lower than that of the disc, this also indicates that 
the continued accretion of gas is driving the change in the 
disc’s angular momentum direction. At tiookback < 1-3 Gyr, 
the central disc decouples and relaxes into the plane of the 
stellar ‘disc’, but the gas at larger radii remains misaligned 
at a constant angle of about 55°. The gas at large radii can¬ 
not be driving the late time relaxation of the inner disc when 
the angle is less than 50°. It is possible that accreting gas, 
which has a larger misalignment angle, slows down the re¬ 
laxation, but it will not be able to speed it up. 

The decoupling of the gas disc and the gas halo hap¬ 
pens after the gas mass has stopped growing, the gas accre¬ 
tion rate has decreased, and the gas metallicity has started 
to increase again, as described in Section 13.11 This is con¬ 
sistent with the picture that in the beginning of the evo¬ 
lution of the misaligned gas disc, the angular momentum 
brought in through gas accretion is high compared to the 
angular momentum change induced by stellar torques. In 
the later stages of the evolution, the gas accretion rate is 
lower and the mass of the disc is higher, resulting in gravi¬ 
tational torques dominating the change of the disc’s angular 
momentum vector. We checked that the total angular mo¬ 
mentum brought in by gas accretion is sufficient to explain 
the angular momentum change of the gas disc in the first 
1.5 Gyr after the merger, but that it is not afterwards. 

In the absence of gas accretion, this gas disc would re¬ 
lax to a misalignment angle of 180°, since it starts out at an 
angle larger than 90°. The fact that this does not happen 
again shows the important role of continued accretion from 
the gas around the central disc, which drives the misalign¬ 
ment angle to less than 90°. 

As explained in the introduction (Section [T]), one would 
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Figure 6. The misalignment angle between the gas and stellar 
angular momentum vectors as a function of lookback time as in 
Figure[5l but dividing the gas disc into 5 shells of 1 kpc thickness. 
The shells at larger radii, R, become misaligned first around the 
time of the merger, while the inner kpc misaligns last. For the 
first part of the evolution (tiookback ~ 3 — 2 Gyr), the shells have 
similar misalignment angles (and the gas disc is not warped). At 
^lookback iS 1-9 Gyr, the smaller the radius of the shell, the faster 
it relaxes into the plane of the stellar ‘disc’, resulting in a warp. 


naively expect relaxation into the plane of the stellar ‘disc’ 
to take only a few dynamical times. However, we find that 
substantially misaligned gas disc is present for about 2 Gyr, 
which is much longer than a dynamical time, even at a ra¬ 
dius of 5 kpc (see Figure [7| below). The persistence of the 
misalignment is mostly driven by continued misaligned ac¬ 
cretion, at least down to tiookback ~ 1-3 Gyr. 

The dark matter angular momentum in the central 
20 kpc is misaligned with that of the stars by about 25° 
during most of the past 5 Gyr. The misalignment with the 
gaseous angular momentum vector is generally larger. There 
is no smooth trend as seen in Figure[5]and we emphasize that 
at these radii the stellar component dominates the potential 
and thus the torques on the gas. The dark matter at larger 
radii also exhibits large misalignment angles with both the 
gaseous and stellar angular momentum in the galaxy centre. 
We are therefore confident that the dark matter is not the 
main driver of the gas-star realignment. 

Figure |6] shows the same misalignment angle as Fig¬ 
ure [S] with the difference that we subdivide the gas disc 
into spherical shells of 1 kpc thickness around the galaxy 
centre. All gas shells become misaligned with the stars 
around the time of the merger, but the shells at larger radii 
misalign first. At tiookback < 1-9 Gyr, the smaller the ra¬ 
dius of the shell, the faster it aligns its angular momen¬ 
tum vector with that of the stars. For example, the in¬ 
ner shell (light-blue curve) drops below 20° misalignment 
at tiookback = 1-3 Gyr, while the outer shell (yellow curve) 
does so only at fiookback = 0.6 Gyr. This causes the gas disc 
to develop a warp, as seen in the fourth and fifth panel of 
Figure[T] The warp is initially most pronounced in the centre 
and at late times in the outskirts. 

Especially the inner shell shows very non-monotonic 
behaviour, with the misalignment angle decreasing and in¬ 
creasing sharply. This indicates a strong role for gas accre¬ 



Figure 7. The relaxation time, determined in two different ways, 
of gas shells divided by the dynamical time at the middle of the 
shell as a function of the radius of the shell. The solid, black curve 
shows the total time the individual shells are misaligned compared 
to their dynamical time, i.e. starting at tiookback = 3.1 Gyr (when 
all gas shells are strongly misaligned with the stars) and ending 
when the misalignment of an individual shell drops below 20°. 
The relaxation time defined in this way does not scale linearly 
with tdyn- The dashed, red curve shows the time it takes each 
shell to relax from 50° to 20°. 50° is the misalignment angle of 
the gas outside the disc (see Figure [5ll, so below this value the 
relaxation cannot be driven by gas accretion. When defined in this 
way, trelax is about 6 tdym expected. This shows that initially 
the gas disc’s evolution is driven by external accretion, whereas 
at late times gravitational torques dominate. 

tion and feedback-driven outflows, since torquing with the 
stellar component will only decrease the misalignment an¬ 
gle. Another possibility is that the disc torques with satel¬ 
lites that pass at a short distance, but this would disrupt 
the outer disc more strongly than the inner disc. 

We refer to the time it takes the gas to realign its angu¬ 
lar momentum direction with that of the stars as the relax¬ 
ation time, treiax- We investigate the relaxation of five 1 kpc 
thick gas shells around the galaxy centre. From Figure [6] it 
is clear that at tiookback = 3.1 Gyr all gas shells have simi¬ 
larly large misalignment angles, of about 120°. We calculate 
Ireiax from this time until the misalignment angle becomes 
less than 20°, when we consider the disc to be sufficiently 
relaxed. The solid, black curve in Figure [7| shows this re¬ 
laxation time (which varies from 1.8 to 2.5 Gyr) divided by 
the time-averaged dynamical time at half the shell radius, 
tdyn = 27ri?i/2/r’circ (wMch varies from 0.01 to 0.12 Gyr), 
where R 1/2 is the radius at the middle of the shell and Ucirc 
is the circular velocity at Ri/ 2 - In the centre the misalign¬ 
ment is visible for over a hundred dynamical times and for 
several tens of dynamical times at larger radii. The relax¬ 
ation time defined in this way does not scale linearly with 

tdyn. 

If only gravitational torques were Important for the re¬ 
laxation of the gas disc, one would expect a constant scaling 
with the dynamical time, which is clearly not what we find 
in Figure [6l However, as we showed above, in this cosmo¬ 
logical simulation, it is not just the torques with the stars 
that are driving the evolution of the gas disc, but also the 
continued gas accretion from its surroundings. In Figure [5] 
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we saw that the evolution of the misalignment angle of the 
gas outside the central disc decouples from that of the inner 
disc when it reaches about 50°. It is likely that afterwards 
the evolution of the misalignment is caused mostly by stel¬ 
lar torques (although external accretion could continue to 
slow down the relaxation). We therefore define treiax of the 
individual gas shells in a second way by calculating the time 
it takes the misalignment angle to decrease from 50° to 20°. 
Note that this is a decrease in misalignment angle of only 30° 
and the ‘full’ relaxation time would be at least a factor of 2 
larger. This second relaxation time (which varies from 0.07 
to 0.58 Gyr) divided by the dynamical time is shown as the 
dashed, red curve in Figure^ We hnd that treiax/tdyn ~ 5—8 
and there is no trend with radius, indicating that, in this 
case, fr eiax scales linearly with the dynamical tim e, as ex¬ 
pected (iTohline et al.l[r982l : iLake &: Normanlll98^ . 


4 DISCUSSION AND CONCLUSIONS 

We have shown results from a cosmological zoom-in simu¬ 
lation of an ETG, which experiences a 1:4 (1:1) merger in 
stellar mass (gas mass) at tiookback = 3.4 Gyr {z = 0.3). 
This merger destroys most of the central gas disc. We have 
studied the subsequent formation of a new gas disc, whose 
rotation direction is misaligned with the stellar rotation. We 
have quantihed its evolution and, in particular, its eventual 
realignment with the galaxy’s stellar angular momentum. 

The misaligned disc in our simulation matches the prop¬ 
erties of observed misaligned gas discs. The star-forming gas 
mass is approximately 3x10® Mq at z = 0. This gas is con¬ 
centrated in the inner kpc and extends to about 1/3 of the 
half-mass radius of the galaxy. These pr operties are consis¬ 
tent with those observe d for such ETGs llYoung et al.ll201ll : 
iDavis et al. |[20n|, [ 2 oT 3 ). The merger also produces a kine¬ 
matically decoupled core of young stars in t he galaxy centre, 
as is observed in some of these systems fe.g. lMcDermid et al.l 

l200dl . 

Our results show that a misaligned gas disc can persist 
for several Gyr, much longer than its dynamical time (about 
150 tdyn in the central kpc). This effect is driven by the 
continued accretion of material from minor satellites, which 
dominates the angular momentum evolution of this system 
for the first > 1.5 Gyr after the initial (major) merger. Once 
the accretion rate has dropped sufficiently, gravitational re¬ 
laxation of the gas disc does occur, over a timescale of about 
6 dynamical times. 

Polar disc galaxies, where an outer disc of gas and young 
stars is misaligned by about 90° from a smaller disc of old 
stars, and warps in outer gas discs have been studied in 
previous simulations (e. g. iBrook et al.l l2008l : iRoskar et al.l 
I2OI0I : ISnaith et al.l I2OI2II . These misaligned gas structures 
are found to form and persist for many Gyr through the 
continued accretion of cold gas with misaligned angular mo- 
me ntum, similar to the simulated ETG we discussed here. 

iLagos et al.l (l2015l f assume in their semi-analytic model 
that when the angular momentum of a dark matter halo is 
misaligned from the stellar angular momentum, a misaligned 
gas disc is present. Our results show that indeed gas (and 
dark matter) in the halo can be misaligned from the stars. 
However, the dark matter in our simulation is not aligned 
with the halo gas either. Furthermore, we showed that it is 


not just the gas-star misalignment angle of the accreting gas 
that matters, but also the rate at which mass (and thus an¬ 
gular momentum) is accreted to counteract stellar torques. 
Eventually, after the accretion rate decreases sufficiently, the 
gas disc relaxes into the stellar plane even though the halo 
gas is still misaligned. 

Davis & Bureau (in preparation) use the observed mis¬ 
alignments of gas discs in ETGs to show that a long relax¬ 
ation time is likely required to explain their misalignment 
distributions. In this work they postulate that a continued 
accretion of cold gas could slow relaxation sufficiently. Our 
results are consistent with this scenario and the relaxation 
tim escale we Hnd matc hes theirs closely. 

ISerra et al.l (120141 ') compare their H l measurements of 
local ETGs (on scales of tens of kpc) to cosmological simula¬ 
tions and find qualitative agreement in the cold gas kinemat¬ 
ics, but too few of the simulated ETGs host misaligned gas 
discs. They speculate that a more sophisticated treatment 
of gas physics and feedback is needed in the simulations. It 
is possible that our simulation matches the observed sample 
more closely, since from Figure [5] we can see that the gas is 
misaligned out to at least 20 kpc for at least 4 Gyr, but a 
larger sample of simulated ETGs is necessary to be sure. 

We caution that from a single example we cannot draw 
strong conclusions about the evolution of misaligned gas 
discs in the general population of ETGs. The simulated 
galaxy we have studied in detail may not be representative 
for the behaviour of all ETGs. It sits in a dense group-like 
environment and is therefore not evolving in isolation. It 
has several satellites and it is 3.4 Mpc away from the cen¬ 
tre of another 2 x 10 ^® Mq halo, also containing a group 
of galaxies. Galaxies in the field could have a very differ¬ 
ent accretion history. Studying high-resolution cosmological 
simulations of a large sample of ETGs in a variety of envi¬ 
ronments will be necessary to quantify how rare or common 
the evolution described in this paper is. 

Low stellar surface brightness features of the galaxy 
merger, such as shells and tidal tails, are visible after the 
merger, but disappear before the gas disc aligns with the 
stellar co mponent, consist ent with merger visibility esti¬ 
mates (e.g. lBell et al.ll200dl . Since the gas-star misalignment 
persists for so long, misaligned gas discs are not necessar¬ 
ily associated with observable merger signatures in the stel¬ 
lar photometry, even when the misalignment was originally 
caused by a merger. 

In our simulation, the misaligned gas disc has a sin¬ 
gle kinematic position angle for the first part of its life, i.e. 
the gas shells in Figure [ 6 ] all show similar misalignments for 
tiookback ~ 3.1 — 1.9 Gyr. At later times, gravitational re¬ 
laxation, which acts on a constant multiple of tdyn, causes 
the gas disc to wa rp once accretion has slowed sufficiently. 

I Davis et al.l ll 201 lll did not hnd any evidence for warps in 
their observational sample of such discs, but with the low 
spatial resolution of thei r data it is not cle ar how strong a 
statement can be made. ISerra et al.l (l2014h found warps in 
35 per cent of ETGs, but on larger scales than studied here. 
Future searches for warps in a statistical sample of galax¬ 
ies with high resolution data can help constrain the bal¬ 
ance between gas accretion-driven misalignment (for which 
no warps are expected) and stellar torque-driven relaxation 
(for which we predict warped gas discs). 

In order to complete such an analysis one requires 2D 
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kinematic maps of gaseous and stellar tracers, such as those 
provided by integral field unit (IFU) spectrographs. Existing 
IFU surveys of galaxies are limited in size, but ongoing and 
future surveys, such as the Calar Alto Legacy Integral Field 
spect roscopy Area survey (CALIFA; iGarcfa-Benito et al.l 
I 2 OI 4 I '). the Sydney-A AO Multi-o bject Integra l field spec¬ 
trograph galaxy survey (SAMI; iBrvant et al.l l2015h . and 
the Mapping Nea rby Galaxies at APO survey (MaNGA; 
iBundv et al.ll201^ . will contain hundreds to thousands of 
objects. These surveys will be able to make exquisite mea¬ 
surements of the kinematic misalignment distribution, and 
how it changes in different objects (e.g. with galaxy mass, en¬ 
vironment etc.). Under the assumption that (major) mergers 
are responsible for creating misaligned gas discs and that we 
know how long the misalignments persist, these observations 
will become a powerful way to constrain the (major) merger 
rate, and its variation, in the nearby universe. If indeed con¬ 
tinued gas accreting significantly affects the timescales over 
which misaligned discs remain observable, these surveys may 
teach us about cosmological gas accretion onto the ISM as 
well. Theoretical studies of the factors that affect the mis¬ 
alignment and realignment of gas discs in galaxies, using 
a large statistical sample of simulated ETGs, will be vital 
to understand the creation and persistence of misaligned gas 
discs and enable correct interpretation of these observations. 
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